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Results are presented from the ATLAS collaboration from the 2010 LHC heavy ion run, during 
which nearly 10 inverse microbarns of luminosity were delivered. Soft physics results include charged 
particle multiplicities and collective flow. The charged particle multiplicity, which tracks initial 
state entropy production, increases by more than a factor of two relative to the top RHIC energy, 
with a centrality dependence very similar to that previously measured at RHIC. Measurements of 
elliptic flow out to large transverse momentum also show similar results to what was measured at 
RHIC. Extensions of these measurements to higher harmonics have also been made, and can be 
used to explain structures in the two-particle correlation functions that had long been attributed 
to jet-medium interactions. New hard probe measurements include single muons, jets and high pr 
hadrons. Single muons at high momentum are used to extract the yield of bosons and are 
found to be consistent within statistical uncertainties with binary collision scaling. Conversely, jets 
are found to be suppressed in central events by a factor of two relative to peripheral events, with 
no significant dependence on the jet energy. Fragmentation functions are also found to be the same 
in central and peripheral events. Updated asymmetry results confirm previous measurements with 
increased statistics, and multiple jet cone sizes are presented. Finally, charged hadron spectra have 
been measured out to 30 GeV, and their centrality dependence relative to peripheral events is similar 
to that found for jets. 



I. HEAVY IONS WITH THE ATLAS DETECTOR 

The ATLAS detector at the LHC [J, shown in Fig. [l] was designed primarily for precise measurements in 
the highest-energy proton-proton collisions, particularly to discover new high-mass particle states. However, 
it is also a very capable detector for the measurement of heavy ion collisions at the highest energies achieved 
to date, already a factor of 14 higher than available at the Relativistic Heavy Ion Collider. The ATLAS 
inner detector is immersed in a 2 Tesla solenoid magnetic field and provides precise reconstruction of particle 
trajectories typically with three pixel layers and four double sided silicon strip detectors out to |7^| = 2.5, as 
well as a transition radiation tracker, out to |7^| = 2. The longitudinally segmented ATLAS calorimeter provides 
electromagnetic and hadronic measurements out to |7^| = 4.9, with particularly high spatial precision in the r] 
direction. Finally, the ATLAS Muon Spectrometer is located outside the calorimeters (which range out most of 
the hadronic backgrounds) and measures muon tracks out to |7^| < 2.7. 

The LHC provided lead-lead collisions to the three large experiments in November and December 2010, and 
ATLAS accumulated almost ten inverse microbarns of luminosity, as shown in the right panel of Fig. [T] The 
analyses shown here typically use up to eight inverse microbarns, during which the main solenoid was activated. 
A smaller dataset with field off was used for multiplicity analyses. 



Muon Detectors Tile Calorimeter Liquid Argon Calorimeter 




FIG. 1: (left) Schematic diagram of the ATLAS detector, showing the three main subsystems: the inner detector 
{\r]\ < 2.5), the calorimeter (|?7| < 4.9) and the muon spectrometer (|?7| < 2.7). (right) Integrated luminosity taken by 
ATLAS in the 2010 heavy ion run vs. time. 
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FIG. 2: Measured distribution of FCal with percentile bins indicated. 



II. GLOBAL OBSERVABLES 



Global observables, such as the total multiplicity and elliptic flow, address crucial issues related to the basic 
properties of the hot, dense medium. The charged particle multiplicity is important both from a thermodynamic 
standpoint, as it is proportional to the initial state entropy, as well as from the perspective of energy loss calcu- 
lations which require an estimate of the initial gluon density. In ATLAS, the multiplicity has been measured [2] 
with the solenoid field turned off to mitigate the losses from the bending of low momentum particles. The pixel 
layers within |?7| < 2 are used to minimize the non-primary backgrounds by only using hits consistent with a 
minimum-ionizing particle emanating from the primary vertex. The pixel hits are assembled both into "track- 
lets" (two hits, both consistent with the measured primary vertex) and full tracks (three points, reconstructed 
with the main ATLAS tracking software). In all, three methods are used, each with quite different systematic 
uncertainties, providing accurate measurements up to the highest particle densities and uncertainties of only 
a few percent. The collision centrality is estimated using the total transverse energy measured in the ATLAS 
forward calorimeter (FCal, covering 3.2 < |7^| < 4.9) and bins are selected using an estimated sampling fraction 
of / = 100 ±2% in the preliminary results shown at Quark Matter 2011, and a revised final value of / = 98 ±2% 
in recently-submitted papers [H [3] . The final centrality selections [3] are shown in Fig. [5] 

The final multiplicity results are shown in Fig. [sj The charged-particle yield, averaged over I77I < 0.5, 
dNch/df]/ (Npartf^) for the 6% most central events, agrees well (consistent well within the stated systematic 
uncertainties) with the previously-available ALICE and CMS data [4 . There is also good agreement for all 
centrality intervals (parametrized by Npart)^ even for the most peripheral events (as shown in detail in the 
inset). After scaling the RHIC data by the observed factor of 2.15, the centrality dependence at the LHC is also 
seen to be quite similar to that observed at a much lower energy. A similar observation was made at RHIC, 
by comparing data from -s/snn = 200 GeV to that from 19.6 GeV and finding a similar centrality dependence 
there as well. Given the dramatic increase in the rates of hard processes at higher energies, this seems to suggest 
that the centrality dependence of the charged particle density is primarily geometric in its origin. 

The collective response of the system to the initial geometric configuration of the nucleons in the colliding 
nuclei is typically characterized by the "elliptic flow", or a cos(2[0 — ^2]) modulation in the azimuthal angular 
distributions relative to the measured "event plane" angle ^2- In ATLAS, the event plane is measured in the 
FCal. Elliptic flow is measured using tracks in the inner detector out to I77I = 2.5. The effects of "non flow" are 
minimized by maximizing the gap between the track and the event plane measurement by using the event plane 
from the opposite FCal hemisphere to the measured track [3 . The half-amplitude of the elliptic flow harmonic 
('^'2)5 as a function of transverse momentum and pseudorapidity, has been measured out to = 20 GeV and 
1 77 1 < 2.5 for eight centrality bins. There is no substantial pseudorapidity dependence out to the maximum 
measured 7^, but there is a strong dependence with a pronounced rise in V2 out to = 3 GeV, a slower 
decrease out to px = 8 GeV and then a very weak dependence out to the highest measured pt- The top panels 
of Fig. [4] shows the final ATLAS data on V2 vs. pt in eight centrality bins. The bottom panel shows the ATLAS 
data from the 40-50% interval compared with already-published ALICE [5^ data as well as lower energy data 
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FIG. 3: (top) Energy dependence of the charged-particle multiplicity density averaged over |?7| < 0.5 for the 0-6% most 
central events compared with heavy ion and proton- (anti) proton data [2 . (bottom) The final centrality dependence of 
the multiplicity per participant pair 2 , compared with ALICE and CMS data and an average of the RHIC data. 



from PHENIX [6^ and STAR [7, all from the 40-50% centrality bin (except STAR data from 40-60%). One 
finds that all of the data are quite similar, even at high within the large statistical errors of the PHENIX tt^ 
data. 

It is not obvious if this apparent scaling behavior, especially at high px, is consistent with expectations from 
jet energy loss. A first look at at high px from energy loss calculations has been made by by Horowitz and 
Gyulassy [8 and is shown in Fig. [s] These calculations are found to describe high px hadron suppression data at 
RHIC quite well. An extrapolation to LHC energies is performed by scaling the initial gluon density by the ratios 
of the charged-particle multiplicity near r] = 0. However, it is found that the hadron suppression thus predicted 
at the LHC is a factor of two smaller than is observed in the recent LHC data (i.e. the measured Raa is a factor 
of two higher than predicted). The predicted values of V2 reflect the differential between energy loss of jets 
emitted "in plane" (i.e. which travel a shorter path length through the medium, and thus are less suppressed) 
and "out of plane" (which are conversely more suppressed). When compared to the preliminary ATLAS data, 
it is found that the data agree surprisingly well with the predictions for > 10 GeV and centrality above 
30%. The more central data (10-20%) show a sigificant divergence from the predictions. However, this may be 
explained in part [9 by the lack of fluctuations in the initial nucleon configurations used in the calculations [10 . 

After the realization that initial state fluctuations can induce higher Fourier modes in the overlap region of 
the two nuclei [11 , elliptic flow has now been joined by the study of higher-order harmonic flow coefficients [12]. 
These are measured each in their own event plane. The resolution has been evaluated by a variety of subdetector 
combinations and it is found that ATLAS can resolve modes up to n = 6 in the most central 40% of events. These 
are shown as a function of pt for six centrality bins in the left panel of Fig.[6j These figures show the rapid change 
of V2 with centrality, expected from previous measurements and attributed to the changing overlap geometry 
as the impact parameter is increased. They also show only a modest change in the higher order coefficients as 
the centrality changes, suggesting that they are mainly sensitive to fluctuations of nucleon positions relative to 
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FIG. 4: (top) Elliptic flow as a function of transverse momentum integrated over |?7| < 1 for eight centrality intervals 
(bottom) Elliptic flow as a function of transverse momentum for the 40-50% centrality interval (^ and references therein), 
compared with data from RHIC and the LHC. 



the overall elliptical shape given by V2. Hydrodynamic calculations suggest a strong dependence of on n, 
particularly when viscous effects are taken into consideration. Thus, it is expected that the information from 
these higher coefficients can distinguish between competing physics scenarios attempting to explain the initial 
state and will provide new information on the viscosity to entropy ratio. 

A longstanding puzzle at RHIC is the observation of unusual structures (e.g. the Mach Cone [13 ) in the two 
particle correlation function, particularly on the "away side" relative to a high particle. Both STAR [14] 
and PHENIX [15 observed a strongly- modified away-side peak which even developed a "dip" at = tt in the 
most central events. The "near side" is also found to show an "ridge" -like enhancement at ~ at large 
values of A77 [16 . However, both of these phenomena have been argued to arise most likely from the presence 
of higher flow coefficients beyond elliptic flow [TT]. In ATLAS, the correlation function has been measured with 
large separation in r] to suppress the contribution from jets and a discrete Fourier transform (DFT) is used to 
extract coefficients out to n = 6 [3]. These are found to agree quite well with the Vn extracted using event 
plane approaches. More interestingly, when the event plane coefficients are used to construct a comparable 
two particle correlation function dN/d(j) oc cos(n(/)), one finds a striking agreement with the correlations 

measured directly, as shown in Fig. |6] This suggests that what were formerly thought to be indications of jet- 
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FIG. 5: Comparisons of ATLAS data of V2 at high pt compared with differential energy loss. 



medium interactions may well result simply from the presence of higher-order flow harmonics, arising through 
a combination of nucleon position fluctuations and viscous effects. 



III. HARD PROBES 



High pt signals produced by perturbative processes in the initial collisions of the two nuclei (also known as 
"hard probes") provide a means to probe the conditions of the collision at very early times. This is typically 
done by means of "suppression" observables which indicate deviations from the expected scaling of yields by 
the number of binary collisions, estimated using Glauber modeling (e.g. the ATLAS J/^l) measurement [E]). 
To confirm the assumed scaling, ATLAS has used the spectrum of single muons at large px to extract the yield 
of bosons, by a simultaneous fit of a template trained on simulated as well as a parameterization of 
muons from heavy flavor quark decays, as shown in the left panel of Fig. [71 A sample of approximately 400 W 
bosons has been extracted. The binary scaling of the measured yields is studied using the variable i^pc, defined 
as the ratio of yields measured in different centrality classes to the yield measured in the 10% most central 
events, with all yields scaled by the corresponding number of binary nucleon-nucleon collisions. The right panel 
of Fig. [rlshows Ry>c as a function of centrality [18 . Using a fit to a constant value, giving (Rpc) = 0.99±0.10 
with a = 3.02 for 3 degrees of freedom, a significant consistency with binary scaling is observed. 

ATLAS has already published striking results on dijet asymmetries [19 . New results are reported here on 
the measurement of inclusive jet yields as a function of jet E^, in order to test models of QCD energy loss 
more directly [20 . Jets are reconstructed using the anti-/ct algorithm with the jet size set to both R=0.2 and 
R=0.4, based on "towers" composed of calorimeter cells integrated over regions of size Arj x = 0.1 x 0.1. 
The ambient background is removed at the cell level by excluding regions near jets and calculating the mean 
energy, as well as any azimuthal modulation, in strips of width Ar] = 0.1. An iterative procedure is applied to 
remove any residual effect of the jets on the background subtraction. The final jets are corrected for the energy 
scale and resolution based on PYTHIA jets embedded into HIJING, with conservative systematic uncertainties 
assigned to account for the small differences between the fluctuations seen in data and simulation. Jets are 
restricted to \r]\ < 2.8, to stay within the main barrel and endcap regions of the calorimeter. 
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FIG. 6: (left) Higher order Fourier coefficients {v2 — vq) as a function of transverse momentum for six centrality bins, 
using the Event Plane (EP) method. It is observed that while V2 varies strongly with centrality, V2, — vq are nearly 
invariant, (right) Measured two-particle correlation function for 2 < < 3 GeV and 77 < 2.5 in the 1% most central 
collisions (open circles) . The reconstructed correlation function using the coefficients from the event plane measurement 
(and vi from the two-particle measurement) is also shown and agrees quite well with the directly measured correlation 
function, using both full FCal (solid line) and FCalp(N) (dashed line) methods. 




p!^[GeV] 1 -centrality 

FIG. 7: (left) Single muon spectrum measured for the 0-10% most central events, with the templates for W bosons and 
heavy flavor (indicated as "Background") also shown to illustrate the yield extraction procedure, (right) Rpc for W 
bosons as a function of centrality, showing consistency with binary collision scaling. The dotted line is a fit to a constant. 



The left panel of Figure [s] shows Rqv for R=0.4 jets, (defined similarly to i^pc, but relative to the 60-80% 
centrality interval) as a function of centrality for jets in three fixed bins (£^t = 50 — 75 GeV, 75 — 100 
GeV and 100 — 125 GeV). A smooth evolution of Rqv is observed going from the 50-60% central events on the 
right to the 0-10% central events on the left, where a suppression of roughly a factor of two is observed for 
all bins and for both R=0.4 (left panel) as well as R=0.2 (shown the right panel). The invariance with R 
is suprising given the general tendency of radiative calculations to predict a modified fragmentation function 
for quenched jets, with substantial transverse broadening. However, direct measurements of the transverse 
and longitudinal fragmentation functions for tracks with > 2 GeV, shown in the left and right panels 
respectively of Fig. [9j confirm that no substantial modification of the fragmentation function can be observed 
when comparing peripheral and central events. In other words, the increased suppression of the jet rates is not 
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FIG. 8: (left) Rcp for R=0.4 jets, as a function of centrality, with 60-80% used as the peripheral sample. In this 
convention, the most central events are on the left and the most peripheral are on the right, (left) Similar as the other 
figure, but for R — 0.2. 
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FIG. 9: (left) Transverse fragmentation function for R=0.2 anti-Zct jets, comparing central 0-10% events with peripheral 
40-80%. (right) the similar comparison for the longitudinal fragmentation function. 



accompanied by any evident modification of the jets themselves. 

The original dijet asymmetry results provided the first direct evidence that jets lost energy through the hot, 
dense medium. The preliminary updated asymmetry results provide several major advances over the earlier 
ones. The full integrated luminosity (Lint = 7 fih~^ is now used, which increases the jet statistics by more than 
a factor of five. For the larger R=0.4 jets, the background subtraction now accounts for the effect of elliptic 
flow, using the measured event plane from the forward calorimeter and the cos(2(/)) modulation measured in 
the different calorimeter layers. Finally, a residual effect of the background subtraction on the reconstructed 
jet energy is removed using an iterative procedure that performs a second background subtraction excluding all 
reconstructed jets with > 50 GeV. This has an approximately 10% effect on jet energies. 

The new results all use a leading jet with Et > 100 GeV and a subleading jet with Et > 25 GeV, with both 
jets within |7^| < 2.8 and > 7r/2. The cut in pseudorapidity is to exclude the ATLAS forward calorimeter, 
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in which the centrahty and reaction plane measurements are made, and the cut in azimuth is to require the jets 
emerge in opposite hemispheres. The top panel of Fig. [To] shows the asymmetry distribution for reconstructed 
R = 0.4 jets, with the elliptic flow subtraction performed, in six centrality intervals. It is observed that the jets 
measured in the most peripheral bins resemble both fully simulated events with PYTHIA dijets embedded into 
HIJING background events, as well as proton-proton data at ^/s = 7 TeV. As the centrality is increased towards 
events with greater nuclear overlap, the distribution measured in lead-lead data becomes broader, and becomes 
essentially flat for the 0-10% most central events. It should be noticed that the slight "peak" structure seen in 
the original results [19 is no longer visible with higher statistics and better control over the background and 
jet energy. The middle panel of Fig. [lO] shows the same results, but for R=0.2 jets. The identical trend is seen 
for these much smaller jets, which have a jet area a factor of four smaller relative to the R=0.4 jets and no flow 
subtraction. The fact that the asymmetry distribution similarly widens for the narrow jets, which are much 
less sensitive to residual background fluctuations, provides a counterargument to the one made in Ref. [21 that 
the asymmetry distribution could be explained by a poor understanding of the underlying event fluctuations. 
The bottom panel of the figure shows the distribution of Acj) between leading and subleading R=0.4 jets. In all 
centrality intervals, the emission remains primarily back-to-back. However, there is a low-lying but significant 
tail at A(j) < 37r/4, which can presumaby be explained in part by the presence of fake jets. 

To study this phenomenon further, charged particle yields have been measured out to |7^| = 2.5 and = 30 
GeV [22 . Tracks are reconstructed in the full inner detector (ID), with stringent requirements on the number 
of hits per track in order to suppress fake tracks. With these requirements, efficiencies are typically 70-80% 
while fake rates near r] = are less than a percent, although both efficiency and purity are degraded somewhat 
going to forward rapidity. While tracks are measured in ATLAS at transverse momenta well past 30 GeV, the 
measured track impact parameter errors (which are also used to reject fake tracks) are not the same in the data 
and in the present ATLAS simulations. 

The top three panels of Fig. [TT] show the corrected spectra out to 30 GeV, in three inclusive pseudorapidity 
intervals (|7^| < 0.35, |7^| < 0.8 and |7^| < 2.5). The statistical errors are shown by the error bars (which are 
typically smaller than the graph markers) and total systematic uncertainties are shown by the light grey bands. 
Dividing the spectra by the yields measured in the 60-80% interval as well as the ratios of the number of 
binary collisions gives the charged-hadron Rcp values shown in the bottom three panels of Figure [TT] for three 
centrality intervals relative to the 60-80% centrality interval. The trend for the 0-5% interval is similar to that 
seen in charged particle Raa measured by ALICE [23j . There is a peak at 2-3 GeV, a minimum around 7 GeV 
and then a steady rise going out to high px which has not yet leveled-off even at px = 30 GeV. The value of 
Rcp averaged over px from 20 to 30 GeV, shown in Fig. [T2j is found to decrease smoothly from a value near 
unity for 50-60% to approximately 0.4 for the 0-10% most central collisions. This trend is both quantitatively 
and quantitatively similar to the one previously discussed for fully reconstructed jets. 



IV. CONCLUSIONS 



Results are presented from the ATLAS detector using over 9 /ib~^ of lead-lead collisions from the 2010 LHC 
heavy ion run. The charged particle multiplicity increases by a factor of more than two relative to the top RHIC 
energy, with a centrality dependence very similar to that already measured at RHIC. Measurements of elliptic 
flow out to large transverse momentum also show similar results to what was measured at RHIC. Measurements 
of higher harmonics have also been made, which are able to explain structures in the two-particle correlation 
(e.g. the "ridge" and "cone" phenomena) previously attributed to jet-medium interactions. Single muons at 
high momentum are used to extract the yields of bosons as a function of centrality, which are found to be 
consistent with binary collision scaling. Conversely, jets are found to be suppressed in central events by a factor 
of two relative to peripheral events, with no signiflcant dependence on the jet energy. Fragmentation functions 
are also found to be the same in central and peripheral events. Updated asymmetry results are presented for two 
different jet radii and improved background subtraction, conflrming the flrst ATLAS results. Finally, charged 
hadrons have been measured out to 30 GeV, and the measured values of Rcp are found to have a centrality 
dependence (relative to peripheral events) similar to that found for jets. 
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FIG. 10: (top) Asymmetry distribution as a function of centrality for R = 0.4 jets (middle) 
a function of centrality for R = 0.2 jets (bottom) Atj) distribution for R = 0.4 jets. 
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FIG. 11: (top) Invariant yields of charged particles for four centrality intervals and three inclusive pseudorapidity 
intervals, (bottom) Charged particle Rcp measured relative to the 60-80% bin, for all tracks measured within |?7| < 0.35, 
\r]\ < 0.8 and |7y| < 2.5. 
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FIG. 12: Rcp for pT = 20 — 30 GeV as a function of centrality, all relative to the 60-80% centrality interval. 



